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Abstract 

The generalization of the Maupertuis principle to second-order Varia- 
tional Calculus is performed. The stability of the solutions of a natural 
dynamical system is thus analyzed via the extension of the Theorem of 
Jacobi. It is shown that the Morse Theory of the trajectories in the 
dynamical system is identical to the Morse Theory of geodesies in the 
Jacobi metric, even though the second-variation functionals around the 
action and the Jacobi length do not coincide. As a representative exam- 
ple, we apply this result to the study of the separatrix solutions of the 
Garnier System. 
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1 Introduction 



The Jacobi version of the Maupertuis principle establishes that the dynamics 
in a natural system can be viewed as geodesic motion in an associated Rie- 
mannian manifold: If 5 = J(T — U)dt is the (natural) action functional of 
a system defined in a Riemannian manifold, (M,g), the critical trajectories 
of S with energy E = T + U coincide with the extremals (geodesies) of the 
length functional = f ds defined in {M,h), where h is the Jacobi metric 
h = 2{E - U)g. 

In this work we analyze whether or not this idea works at the level of 
second-order variational calculus. Namely, can the Hessian operator of a nat- 
ural dynamical system be viewed as the geodesic deviation operator of the 
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associated Riemannian manifold? Early on, it was pointed out that the an- 
swer is negative Q and consequently the stability criterion provided by the 
geodesic deviation operator ("geometric stability criterion") does not coin- 
cide with the usual criterion based on analysis of the spectrum of the Hessian 
operator ("dynamical stability criterion"). We shall show that even in this 
situation the Morse theory associated with the dynamical problem is the same 
as the Morse theory built on the geodesies of the Jacobi metric. The crux of 
the matter is that the Jacobi weak stability criterion provided by the number 
of conjugate points crossed by a trajectory is the same from both points of 
view. 

There are several topics where the geometric criterion has been used, 
e.g. study of the chaotic behaviour of Hamiltonian systems Q or the non- 
integrability of dynamical systems [Q. The Jacobi-Maupertuis principle has 
been extended to the case of Lorentzian manifolds in [12|. 



2 The Maupertuis-Jacobi Principle 
2.1 Preliminaries and Notation 

Let M = (M, g) be a Riemannian manifold. Geodesies in M are extremals 
of the free-action So or length functionals L for any differentiable curve 7 : 

Soil] = r \m)fdt; m = r \m\\dt . (i) 

The critical paths are the solutions of the Euler-Lagrange equations: 
J5o = O^V^7 = 0; 5L = ^ V^7 = A(t)7, = ( ^ ) 

(2) 

As a consequence of the invariance of the length functional under reparametriza- 
tions, variational calculus on L leads to equations where the geodesies are 
parametrized by an arbitrary parameter t (often called pre-geodesics). 

To decide whether a critical path is a local minimum we focus on the 
second variation functionals: 

{^V,V)ds, 6^L= /AV^,V^\ds 
AV = -VyVy V - R{i,V)i = - K^'{V) ■ 
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Here, V denotes a proper variation {V{si) = ^^(52) = 0, see for instance Q 
for details) and V-^ is the component of V orthogonal to the geodesic. R and 
K are respectively the curvature and sectional curvature tensors. Thus, A is 
the geodesic deviation operator. 

In a natural dynamical system where the configuration space is the Rie- 
mannian manifold {M,g), the dynamics is governed by the action: 

%]= j[*'Qll7lP-f/(7W)) , (3) 

and the Euler-Lagrange-Newton equations are: 

5S = 0^^ = -gradf/ . (4) 
dt ^ ^ ^ 

The local stability of a given extremal under proper variations is determined 
by the Hessian or second-variation functional: 

S'^S = r dt {{AV, V) - H{U){V, V))= r dt {AV - Vygradf/, V) 
Jti Jti 

(5) 

and we shall denote the differential operator in the quadratic form ^ as: 
AV = AV- VygradC/. 



2.2 Geodesies in the Jacobi Metric 

Let X* be a system of local coordinates in M, and let us set ds'^ = dij^-x"^ ® 
dx^ = dijdx^dx^ . The Jacobi metric h is defined as a conformal transfor- 
mation of g: ds\ = hijdx''dx^ , hij = 2{ii — U)gij, where ii is a constant 
and the Riemannian character of h restricts the admissible values for ii by 
means of the inequality U < ii. 

Given any two vector fields X,Y €z T{TM), we shall use the following 

notation: h{X,Y) = {X,Y)-^ , \\X\\-^ = ^ {X^X^ , and we shall write s for Sh 
for simplicity: ds = dsh- 

Theorem of Jacobi. The extremal trajectories of the variational problem 
associated with the functional with energy E = ii, are the geodesies of the 
manifold (M, h), where h is the Jacobi metric: h = 2{ii — U) g. 

The proof of this theorem can be found in several References (see Q for 
instance). From an analytic point of view, the theorem simply establishes 
that the Newton equations @) arising from the action S are tantamount to 
the geodesic equations in (M, /i):Vy7' = 0, 7' = The equivalence is 
shown by performing both a conformal transformation, h = 2{ii — U)g, and 
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a reparametrization from the dynamical time t to the arc-length parameter s 
in (M,/i), 

^ = 2V(ii-^(7(^)))r = (n-^(7(.))) , (6) 

over the solutions (note that in the domain U < ii the reparametrization is 
well defined). 

3 "Geometric Stability" versus "Dynamical Stabil- 
ity" 

In the language of Variational Calculus, the Jacobi theorem states that 5S = 
<^ 6Sq = 0, or, for geodesies parametrized with respect to the arc-length, 
6S = <^ 6L^ = 0. In [|^ we showed the following theorems referring to the 
second- variation functionals: 

Theorem 1. Let ■^{t) he an extremal of the S[^] functional, and let [7] 
he the free-action functional for the Jacohi metric associated with 5 [7]. The 
corresponding Hessian functionals satisfy: 

6'Si[j] = 6'S[j] + dt 2 (^7, ^) {F, V) , (7) 

where F = gradln(2(ii — U)). 

Theorem 2. Let 7 he an extremal of the S[y] functional and let L'^l'j] be the 
length functional of the Jacobi metric associated with S['y] . The corresponding 
Hessian functionals satisfy: 

6'lJ[j] = 6'S[^] - / ^ [(V^7, V) - (7, V^V)f (8) 
Jti ^{ii-U) 

Equation (^) shows that geodesies of minimal Jacobi length ^'^[7] are equiva- 
lent to stable solutions of the dynamical system, although the converse is not 
necessarily true. Allowing only orthogonal variations, V = V'^, @ can be 
re-written as: 

The proofs of these two theorems are based on the behaviour of the co- 
variant derivatives and the curvature tensor under reparametrizations and 
conformal transformations of the metric tensor. 
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4 Jacobi Fields and Morse Theory 

A basic tool for the analysis of the stability properties of a given extremal of 
a functional is the study of the Jacobi fields: the elements of the Kernel of 
the second-variation or Hessian operator. Direct computation of the Jacobi 
fields in a given problem is in general a very complicated task: one needs to 
solve the Euler-Lagrange equations, calculate the second- variation operator 
and, finally, determine the Kernel. Fortunately, there exists a short-cut if one 
knows a parametric family of extremals (see, for instance, for details): 

Proposition. Jf 7 = ^{t, a) is a family of extremals of S[y] characterized hy 
the value of the real parameter a, then the vector field ^ is a Jacobi field. 

Knowledge of the Jacobi fields allows us to calculate the conjugate points 
along an extremal 7, and it is thus possible to apply the Jacobi criterion for 
weak minimizers of a functional which counts the number of conjugate points 
crossed by the extremal as a measure of the instability. 

Morse Theory establishes a link between the topology of a manifold and 
the critical-point structure of the functions defined in the manifold [0]. We 
shall deal with Morse Theory in spaces of closed paths with a fixed point in 
Riemannian manifolds. The aim is the application to the separatrix trajecto- 
ries in the Garnier System. A formulation of this theory a la Bott Q is as 
follows: 

Let nM = {7:5^^ M/7(0) = 7(1) = mo} be the loop space in M with 
a fixed base point. Then, some of the topological properties of 0,M are codified 
in the Poincare series Pt (OM") = YlT=o t^, where bk = dim Hk (OM, M) are 
the corresponding Betti numbers. 

For any functional S defined on QM, we define the Morse series as A4t {S) = 
J2n, Pt {Nc) tf'^^-^ where the sum is over the critical manifolds Nc, formed 
either by a single isolated critical path or by a continuous set of critical paths. 
/i(7c), the Morse index, is the dimension of the subspace where the Hessian 
operator of the functional S, A, along 7c is negative definite. 

The Morse inequalities, AdtiS) > Pt (OM), tell us that the topology of 
r^M forces the existence of most of the extremals of the functional S. 

The Morse index can be computed using the Morse index Theorem: 
The Morse index of a critical path 7^ is equal to the number of conjugate 
points to the base point crossed by 7c counted with multiplicity. Note that 
this theorem is equivalent to the Jacobi criterion for weak minimizers. 

From this point of view, knowledge of a one-parametric family of solutions 
7(t; a) of the dynamical problem (P) informs us about the geodesies 7(5; a) in 
{M,h). There is a one-to-one correspondence between the trajectories of the 
dynamical system and the geodesies of the Jacobi metric, which only differ by 
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proper reparametrizations. Therefore, the two kinds of extremals share ah the 
conjugate points, and the structure of the critical points of S and L"^ is the 
same. The Morse series coincide, Ai{S) = Ai{L'^), and the Morse theories 
built from the trajectories in (M, g) and the geodesies in (M, h) are identical. 

5 Separatrix trajectories in the Garnier System 

As an application, we shall analyze the separatrix trajectories in the Garnier 
System The configuration space is (M^,(5ij) and the action functional 
reads: 




Here ai are constants, such that fij ^ aj, Wi ^ j, and the Euler-Lagrange 
equations are: 

^ = 2q,{q-q-l)+ajq^, i = l,2,...,N . (10) 

The system is completely integrable and, in fact, Garnier found all the 
periodic solutions in terms of hyperelliptic functions. In [Q| we computed the 
separatrix solutions lying on the boundary between bounded and unbounded 
motion . We shall focus in the N = 2 case for simplicity (in [l| the very rich 
= 3 case is analyzed thoroughly). 

For any path q{t) = {qi{t), q2{t)) in (R^, (5jj), the Lagrangian C = C{qi, q2, 
91,92) is: 

^ = ^(9? + 9|) + ^(9? + 9i-l)' + y9i . (11) 

With no loss of generality, we have set ai = 0, a2 = cr. We shall only deal 
with the case a < 1, because the separatrix solutions in this regime have a 
richer structure. In fact a > 1 generates a situation where the separatri- 
ces are reduced to a single singular solution. The system is not only com- 
pletely integrable, but Hamilton- Jacobi separable in Jacobi elliptical coordi- 
nates: ql = CT~^(1 — ^i)(l — /i2), 92 = '^~^(^^~/^i)(^^~/^2), where = 1 — cr^. 
With this change, is mapped into the infinite parallelogram determined by 
the inequalities: —00 < ^ui < ct^ < ^2 < 1- 

The choice of the separatrix trajectories is achieved by setting the val- 
ues of the two invariants in involution to zero: ii = i2 = 0. ii denotes the 
energy first-integral and i2 is the second first-integral of the system, function- 
ally independent of ii (see for instance ||To[| for details). This choice confines 
the dynamics to the finite parallelogram P(0) = < /ii < cr^ < /i2 < 1 (in 
Cartesian coordinates it is the area bounded by the ellipse qf + ^ = 1). 
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In elliptical coordinates the Euclidean metric 6ij becomes : gi2 = 
obtained from 



521, 511 = 4i^,-i)(f,%2y 922 = 4(^^l^i)(^7J^2) . The Christoffel symbols are 



pi _ 1 2Ati-(l+o-2) . -pi _ -1 pi _ (A'i-1)(mi-o-^) 

^ 11 - 2(a»i-M2) 2{/ii-l){Mi-<?^)' 12 - 2(Mi-Ai2) ' 22 - 2(m2-1)(M2-S^)(mi-/'2) 

Replacing 1 by 2, is equivalent to and the same happens between Ffg, 
rf^, and r}2, ^22- The rest of the symbols are trivially related to these by 
symmetries or are null. 

The Jacobi metric associated with the system, with ii = 0, is /i = 
fj,2-i_i^ (— (//f — cT^/if) + (/^l — cT^/il))^, ^-iid the Christoffel symbols can be ob- 
tained by direct calculation (||l|). 

5.1 Singular Solutions 

There are two singular solutions to the system, and these live on the border 
of the parallelogram P(0). 

Taking q2 = 0, the equations of the motion can be easily integrated, 
even in Cartesian coordinates, to find the solutions: q{t) = (gi(t),0) = 
(± tanh(t — to); 0). In elliptical coordinates, these solutions lie on two edges of 
the border of P(0): I, fli{t) = (l-tanh^(±t),a-^), fort G (-oo, arctanh(=FfT)] U 
[arctanh(ib(T), oo), and II, flii{t) = (ct^,! — tanh^(=pt)) for t G [arctanh(=F(7), 
arctanh(ib(T)] (taking to = 0). 

The geodesic equations of the Jacobi metric reduced to the ^2 = orbit 
become a single non-trivial differential equation, which is solved by a cubic 
algebraic equation : 

(li + -2 Til = ^ - — = ±s . 12) 

q{-l 3 

(ti, t2) = (— OO, oo) has been reparametrized to [si, S2] = [— |, |] • The explicit 
solution of the cubic (p^) is 

9 

9i(s) = 9i(s) = -cos- + \/3sin- , (13) 

where 9 = arctan ^^"^^'^^ . 

There exists a second type of singular solutions living on the edge /ii = 

0, which also belongs to the border of P(0). This edge becomes the el- 
2 

lipse + ^ = 1 in Cartesian coordinates. The solutions are either q{t) = 
(tanh(ibo"t), cj sech(ib(Tt)) in M^, or /2(t) = (0,1 — cj^ tanh^ crt) in the elliptic 
plane. 
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qi(t) 



qi (s) 




Figure 1: q2 ~ singular solution as a function of the dynamical time t and as a 
geodesic parametrized by the arc- length parameter of the Jacobi metric s. 

2 

The geodesic equations of the Jacobi metric reduced to the q\ + ^ = 1 
orbit also become a single non-trivial differential equation, which is solved by 
a cubic algebraic equation : 

^2.,, . 

s . (14) 



—q-^ =0 ^a[qi 



The explicit solution of the cubic together with the elliptic orbit afford the 
geodesic: 



qi{s) = -i cos f + ^ sin f , ^2(5) = 1^-2 + + cos f + ^/3sin 

2 2 

from si = — (t(1 — ^) to S2 = cr{l — ^). 



29 




- - - - - - 1 








si 








— 1 





Figure 2: Singular solution over the ellipse as a dynamical solution and as a geodesic. 
5.2 General Solution 



The Hamiltonian associated with the Lagrangian (11), in elliptical coordinates 
is of the Stackel type: 
1 



- (-4(^1 - l)(/ii - a'')T:\ - 4(1 - ^2)(/i2 - ^')vr2- 
2(/ii - 112) 

-(/i?-aV) + (Mi-t^V2)) 



(15) 
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and the integration of the Hamilton-Jacobi equation, for separation constants 
equal to zero, [|^, leads to the equation of the orbits, 

sign(7ri) 



a 



sign(7r2) 



/W2 + cr 

and the equation of the temporal dependence, 

o-2sign(7ri) /rj — ^ o-2sign(7r2) 



a 



,2(t+to)a 



(16) 



(17) 



Ai2 + cr 

Here a is a constant that parametrizes the different orbits and to is another 
constant coming from the time-translation invariance. 

At the limit a — > ±oo, equation ([T^) reduces to the singular solution (?2 = 
together with the singular solution living on the "ellipse" . Thus, identifying 
a = —oo with a = oo, we see that all the orbits are parametrized by a periodic 
parameter a : the space of all orbits is the S"^ circle. 




Figure 3: Graphics of the orbits ( p^ in the Cartesian plane and in the elliptical one. 
The graphics were obtained by numerical calculations. 

A straightforward application of the Hamilton-Jacobi procedure to the 
free-Hamiltonian, = ^ (/i^^vf^ -|- /i^^vi"!) , vfj = ^j-, provides all the geodesies 
of the Jacobi metric with ii = 0. 

The equation giving the orbits is exactly (16), and the arc- length depen- 
dence is determined by the expression: 

Ml + 2 n — ^ -,^sienf7r9■)M2 + 2 



s+so = -(-i)-s°(-i)^ii^yr^- (-i)-s-(-2)^i£^yr^ (18) 

5.3 Jacobi fields 

In Q the Jacobi fields are calculated on the trajectories. Here, we closely 
follow the same idea to compute the Jacobi fields on the geodesies, which 
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are the solutions of the equations (J16|)-(pL8D. Thus, we have a two-parametric 
family jl{s; a, sq) of geodesies, and, accordingly, two Jacobi fields: ^ and 
see Section §2.2 . The first Jacobi field is tangent to the geodesic because 
it obeys the invariance under s translations. We focus on the computation 
of the orthogonal Jacobi field: An explicit formula for the geodesic 

fl{s;a,SQ) = (/ii,/U2) in terms of known analytical functions is not available. 



but we implicitly derive the system (1(:)-(|18|) with respect to a to obtain a 
linear system in the components of J: 



-(-i)"i dfii (-i)-2 dfi2 



fii\/l - /Ui((t2 - /ii) da /X2VI - Ai2(o-^ - M2) da 



(19) 



(-irm dfi, ^ (-l)"^/22 dfi2 ^ Q ^20) 
\/l - III da -y/1 - /U2 da 

The Jacobi field that solves (|T9|)-(pO|) is: 



(21) 



J = j(/^i, A2) (^(-i)"iA2 yr^^ + (-i)"^Ai v"!^^) 

where 

-X 2/2i//2(/il -^^)(/i2 -0^2) 

j(m,^2) = 7Z _ w_2 , -2 , - f- , - N-2V "1,02 = 0,1 

(/^l - A^2)(A^i +fJ'2 + ^1/^2 - (/il + ^J'2)(TV 

J is zero both at the starting point D and at the focus F of the ellipse, 
and, henceforth, the focus is a conjugate point of D of multiplicity 1. 



5.4 The Morse series 

Taking M = P2(0); ^IM = C, we have the following critical point structure 
for the geodesies that start and end at the point D. Firstly, the point D itself 
is a possible geodesic with a Morse index 0. Secondly, all members of the 
family determined by the general solution, /z(s; a, sq) cross the focus. Bearing 
in mind that Pt{S^) = (1 + t), the family contributes to the Morse series with 
{l + t)t. Iterating two of the solutions, fl^jl, D is a conjugate point to itself and 
the contribution to the Morse series is (1 + t)t'^. More iterations tell us that 
MtiL-^\c) = l + {l + t)t + {l + t)t^ + ... = ji^. The Morse series is equal to the 
Poincare series of the space of closed geodesies in S'^: Ait = Pt 

In |2| we calculated the Morse series for the = 3 case. The result is: 
Mt{LJ\^)=Pt{C) = j^ = Pt{^S^). 
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